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ABSTRACT 

We present intrinsic extinction curves for 14 AGNs. The AGNs have reddenings, E(B-V), of up to 
0.36 mag. The majority (13 out of 14) of the extinction curves are not steep in the UV. Of the seven 
best determined extinction curves, five have extinction curves that are as flat as the standard Galactic 
curve in the optical and near UV, but flatter in the far UV, and without the A2175 feature. One 
AGN, B3 0754+394, has a steep SMC-like extinction curve, and another, Mrk 304, has an LMC-like 
extinction curve, including a probable A2175 bump. The remaining seven, lower-quality, extinction 
curves have overall shapes that are consistent with an LMC-like shape or a flatter shape. Two have 
possible A2175 features, and one might be identical to the Galactic curve. The flatter curves that 
predominate in our best determined extinction curves are not as flat as the Gaskell et al. (2004) 
extinction curve for radio-loud AGNs. This suggests that the previous radio-loud extinction curve 
might be slightly too flat in the range 4 < 1/A < 6.5 /im~^ because of luminosity-dependent reddening 
biases in the composite spectra, but further investigation is needed. We present a parameterized 
average AGN extinction curve. Observed variations in the continuum properties of the AGNs are 
inconsistent with intrinsic object-to-object variations because observed differences are least in the far 
UV where changes in the accretion disk spectrum should be greatest. We suggest that the steepening 
of AGN spectra around Lyman a is the result of a small amount of SMC-like dust (E(B-V) ^ 0.03). 
We find the largest object-to-object differences in spectral shape to be in the Fe II emission of the 
"small blue bump" . 

Subject headings: galaxies: active — ISMidust, extinction - galaxies:quasars:general — black hole 
physics — accretion: accretion disk 



1. INTRODUCTION 

Knowing the reddening of the continua and emission 
lines of AGNs is vitally important for understanding the 
energy-generation mechanism of AGNs, the physical con- 
ditions of the line-emitting gas, the differences between 
AGNs, and AGN demographics. 

Gaskell et al. (2004) showed that reddening is the ma- 
jor cause of differences in optical to UV spectral energy 
distributions (SEDs), and that the mean extinction in 
AGNs increases with decreasing luminosity. Gaskell et 
al. estimated mean AGN extinction curves from com- 
posite spectra based on samples of AGNs with different 
mean extinctions. These extinctions curves differ from 
the well known extinction curve for dust in the solar 
neighborhood by lacking the A2175 carbon feature. The 
mean extinction curve for the Small Magellanic Cloud 
(SMC) also lacks the A2175 feature, but the mean AGN 
extinction curves found by Gaskell et al. (2004) are con- 
siderably flatter than the SMC curve in the ultra-violet. 
Czerny et al. (2004) also derived a flat extinction curve 
with no A2175 bump from SDSS composite AGN colors. 

The extinction curves derived by Gaskell et al. (2004) 
and Czerny et al. (2004) refer to samples of AGNs. Us- 
ing composite spectra has the advantage of minimizing 
the effects of possible real differences in the spectral- 
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energy distributions, but the composite spectra neces- 
sarily combined quasars of different redshifts. For an 
apparent-magnitude-limited sample the objects at higher 
redshift have higher mean luminosities. The luminosity- 
dependence Gaskell et al. find for the mean extinction 
therefore introduces the possibility of bias in deriving a 
extinction curve for AGNs from composite spectra, since 
the UV part of the extinction curve will tend to come 
from AGNs with lower mean extinction (Willott 2005). 
In this paper we avoid this bias by deriving extinction 
curves for individual AGNs. 

The derivation of the extinction curves is discussed in 
sections 2 and 3, and the curves are presented in section 
4. We consider luminosity dependence in section 5 and 
in section 6 we show that the extinction curves we derive 
cannot be the result of intrinsic variation of the under- 
lying quasar spectral energy distributions or of varying 
amounts of host galaxy light. We discuss our results in 
section 7. 

2. SAMPLE 

To investigate the shape of the extinction curve one 
needs carefully calibrated spectra covering the whole ul- 
traviolet and optical spectral regions. Shang et al. (2005) 
present spectrophotometry of 17 AGNs taken with the 
Far Ultraviolet Spectroscopic Explorer (FUSE) satellite, 
the Hubble Space Telescope (HST), and the 2.1-m tele- 
scope at Kitt Peak National Observatory. The spectra 
cover rest-frame wavelengths of 900 to 9000A and are 
thus well suited for studying reddening. Full observa- 
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Fig. 1. — Spectra of the FUSE/HST sample. Spectra are shown 
in Fx magnitudes normahzed to 1.8 /im~i. Spectra have been 
progressively offset by -1 mag to avoid overlap. Data from Sheng 
et al. (2005). 

tional details can be found in Shang ct al. (2005). 

It is important to note that though the Shang ct al. 
sample has the necessary wavelength coverage for de- 
termining reddenings, it is a very biased sample. The 
objects were specifically chosen with the hope that they 
would be bright in the far-UV (Kriss 2000). This is go- 
ing to bias the sample in favor of low-reddening AGNs. 
Nonetheless, we shall see below that quite a number of 
the AGNs have substantial reddening. 

We show all the spectra used in Fig. 1. We do not con- 
sider the spectral region shortwards of Lyman a. because 
most of the AGNs show absorption lines in this region 
(see Fig. 1 of Shang et al. 2005). All spectra have been 
normalized to the same zero point at 1.8 /xm"^ and cor- 
rected for Galactic reddening. To show the diversity of 
shapes better we show smoothed versions of the spectra 
in Fig. 2 without the offsets. Note that the smoothing 
is purely to reduce confusion for the eye in following in- 
dividual spectra in this figure, and that the smoothed 
spectra were not used in the analysis. 

3. DETERMINING THE EXTINCTION CURVES 

As this paper only addresses the intrinsic reddening 

in the AGNs, we first removed the reddening due to 
dust in the Milky Way. We used the reddening estimates 
of Schlegel, Finkbeiner, & Davis (1998) and a standard 
Milky Way extinction curve as parameterized by Wein- 
gartner & Draine (2001). 

The standard method for determining extinction 
curves is to take the difference in magnitudes between 
a more reddened object and a less reddened object of 
the same intrinsic spectral shape. This is most com- 
monly done for pairs of hot stars of the same spectral 
class. Extinction curves are conventionally given relative 
to the extinction in the V band and normalized so that 
the difference in the extinction in the B band relative to 
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Fig. 2.- - Superimposed spectra of the FUSE/HST/optical sam- 
ple. Spectra arc shown in magnitudes normalized to 1.8 ^im^^. 
The spectra have been smoothed to make it easier for the eye to 
follow individual spectra. The upper of the two thicker lines is 3C 
273 and the lower is B3 0754+394. Data from Sheng et al. (2005). 

the V band is one magnitude (i.e., they are normalized 

to E(B-V) = 1). The scaling factor which the observed 
magnitude differences relative to the V band needs to 
be divided by to give the normalized extinction curve 
is the difference in reddening, E(B-V), between the two 
objects. 

Crenshaw et al. (2001) derived a extinction curve of 
NGC 3227 by comparing it with NGC 4151, and Cren- 
shaw et al. (2002) did the same for Ark 564 by compar- 
ing it with Mrk 493. Since we have many AGNs here, 
and because theory suggests the possibility of real dif- 
ferences between the intrinsic spectra of different AGNs, 
we used an average of several AGNs so that the extinc- 
tion curves were not just reflecting spectral peculiarities 
in one AGN at one epoch. We choose as our templates 
the three AGNs that had the highest relative flux from 
4 to 8 /im~^ compared with their optical fluxes: 3C 273, 
PG 1100+772, and PG 0953+414. Although it was not 
a factor in their selection, these three AGNs happen to 
represent all three main radio classes: 3C 273 is a core- 
dominated source, PG 1100+772 is lobe dominated, and 
PG 0953+414 is radio quiet. 

The extinction curves in Fig. 1 and Fig. 2 of Gaskell et 
al. (2004) were normalized to the B and V bands in the 
standard way, as described above. However, the uncer- 
tainties in normalizing the curves to a relatively small re- 
gion in the optical are amplified considerably by the time 
one gets to the UV, so we adopted a different approach. 
The differences between possible extinction curves, such 
as the standard Galactic curve or the SMC curve, up 
to 3 iJ,m~^ are, for our purposes, very small, so a more 
reliable procedure is to normalize the curves over the 
whole region from the V-band to 3 /im^^. Unfortunately, 
the region around 3 /im^^ is in the middle of the strong 
broad-line region Balmer continuum and Fe II emission 
making up the long-wavelength side of the "small blue 
bump" where there are significant object-to-object vari- 
ations (see Fig. 1). We therefore moved the baseline of 
our normalization to the region of 3.75 to 4.07 /um~^. By 
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TABLE 1 
Estimated Reddenings for AGNs 



wrtPG1 100+772 
wrtPG 0953+414 
wrt 3C 273 

Milky Way minus 2175 
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Fig. 3. — Normalized extinction curves for PG 1351+640 derived 
with respect to the three template blue AGNs. A Milky Way ex- 
tinction curve with the A2175 carbon feature removed is shown for 
comparison. 

4 /ini^^, known, wcU-dctcrmincd extinction curves begin 
to deviate slightly so we experimented with slightly dif- 
ferent normalizations covering the range from SMC ex- 
tinction curves to the Galactic extinction curve for our 
local solar neighborhood. 

In Fig. 3 we show the differences in the separate ex- 
tinction curves of one AGN, PG 1351-1-640, with respect 
to each of the three separate template blue AGNs. This 
shows the reader the range of uncertainty caused by using 
a single template AGN. Note that the differences between 
the curves would be much smaller had we normalized be- 
tween 4.5 and 6 /zm"^. For the rest of the paper, except 
where noted, we will just give extinction curves derived 
with respect to the average of all three template AGNs. 

In our extinction curves we have omitted the regions 
where there are absorption lines. In some spectra there 
are problems with joining up the UV and the optical (see 
Shang et al. 2005), and the signal-to- noise ratio is gener- 
ally worst in these regions. In addition, as can be seen in 
Fig. 3, matches are poorest around the strong emission 
lines. We identified three problems in the region of the 
strong UV emission lines. First, the equivalent widths 
of the strong emission UV lines in AGNs decrease with 
increasing luminosity - the well-known "Baldwin effect" 
(Baldwin 1977, Baldwin et al. 1978) - so to determine 
reddenings from the emission lines one would have to 
match AGNs in luminosity. The second problem is that 
strengths of the cores of the lines also vary with the so- 
called "eigenvector 1" (EVl, Boroson & Green 1992), 
so to determine line reddenings one would also have to 
match AGNs in EVl . The third problem is that the pro- 
files of the UV lines differ from object to object because 
of the blueshifting of the high-ionization lines (Gaskell 
1982). All three of these effects occur to some degree in 
the AGNs we consider. Since we are confining ourselves 
in this paper to the issue of the shape of quasar extinction 
curves, in the majority of cases we have simply deleted 
the parts of our extinction curves around the major UV 
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emission lines. 

4. EXTINCTION CURVES 

The process of normalizing the extinction curve for 
each AGN produces an estimate of E(B-V). We give 
our E(B-V) estimates in Table 1. Column 2 gives E(B- 
V)a2500: the reddening from the normalization as de- 
scribed above. Column 3 gives E(B-V)ai700i the redden- 
ing obtained by comparing the fluxes in the 4 to 8 /Lim~^ 
region with the optical flux, and column 4 gives the av- 
erage of the two estimates. E(B-V)a2500 assumes that 
the Fell emission around A2500 in the small blue bump 
is the same in all AGNs, while E(B-V)ai7oo makes the 
assumption that all the extinction curves are the same 
shape. Our method of deriving extinction curves auto- 
matically makes the average reddening of our template 
AGNs zero. In Table 1 wc have added a small constant 
to reddenings of each set to make the reddenings of the 
bluest template AGN be zero. For the three AGNs with 
very low reddening (see below), E(B-V)a2500 was found 
by comparison with 3C 273 alone. The good agreement 
between E(B-V)a2500 and E(B-V)ai7oo suggests that the 
1-sigma uncertainty in each E(B-V) estimate is ±0.018 
mag, and that the uncertainty in E(B-V)average is ±0.013 
mag. 

Because the quality of derived extinction curves goes 

down as the reddening decreases, we present the extinc- 
tion curves in three groups in order of decreasing redden- 
ing and increasing uncertainty. The groups correspond 
to AGNs with well-determined extinction curves, those 
with less certain extinction curves, and AGNs with little 
reddening and hence ill-determined extinction curves. 

All extinction curves are known to be smooth in the 
UV, and, apart from possible A2175 features, we are 
only interested in this smooth overall shape. We have 
therefore, as discussed above, removed the effects of mis- 
matches in the equivalent widths of the strong emission- 
lines, occasional absorption lines, and, in some cases, the 
poorer signal-to-noise ratio at the junction between the 
optical and HST observations. 

4.1. Well-Determined Extinction Curves 
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Fig. 4. — Extinction curves, derived using the three template 
AGNs, for the six AGNs in the sample with the highest extinction. 
The three smooth curves in each figure are, from top to bottom, 
the SMC curve, an LMC curve, and the standard Galactic curve 
for dust in the solar neighborhood. 

In Fig. 4 we show the extinction curves for the six 
reddest AGNs in the sample. All six lack the A2175 fea- 
ture. Only one AGN, B3 0754+394, has an SMC-Ukc 
curve. Four of the six have curves that arc flatter than 
the Galactic curve in the far UV, and one AGN, Ton 951, 
is consistent with an LMC-like curve, although, again, it 
is flatter in the far UV and without the A2175 feature. 

Quasars with low reddening are most sensitive to differ- 
ences in the spectral shape. These differences are largely 
due to the broad-line region. Our most heavily-reddened 
quasar, PG 1351+640 has significant differences in the 
emission-line strengths and profiles, and also strong ab- 
sorption (e.g., in C IV), but these differences cause rela- 
tively minor differences in the extinction curve (see Fig. 
4). One of the well-known differences bctwccin AGNs is 
the strength of the optical and near UV Fe II emission. 
In the extinction curves one can see the effect of the vary- 
ing Fe II strength especially on the long-wavelength side 
of the "small blue bump". This causes the extinction 
ctu vcs around 3/im~^ to be too high or too low in some 
objects. 

4.2. Less Certain Extinction Curves 

Fig. 5 shows, in order of decreasing quality, the ex- 
tinction curves for the five next most reddened AGNs. 

At lower reddening, the spectral differences mentioned, 
notably the broad Fe II emission which is not easy 
to remove, become more important. The best of the 
curves, that for Mrk 304, is consistent with an LMC-like 
curve including a probable A2175 bump. The other four 
curves are of poorer quality. The extinction curve of PG 
1322+659 is quite close to a Galactic extinction curve 
with a possible A2175 bump, although the low extinction 
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Fig. 5. — Extinction curves for AGNs of lower extinction. As 
in Fig. 4, these arc clerivecl using the three template AGNs. The 
smooth curves again represent SMG, LMC, and Galactic extinction 
curves. 



makes mismatching of emission features a serious prob- 
lem. PG 0947+396 is consistent with cither an LMC-hke 
or a Galactic curve, depending on the normalization. For 
both PG 0947+396 and PG 1322+659 there is probably 
a A2175 bump. The two remaining curves, PG 2349- 
014 and 4C+34.47, are most consistent with LMC-like 
curves, but the flatter shape favored by the majority 
of the AGNs in Fig. 4 is also possible with a different 
normalization. For 4C+34.47, in addition to the uncer- 
tainty in the normalization caused by spectra differences 
around A2500, there is an additional uncertainty caused 
by a gap between the HST spectrum and the Kitt Peak 
spectrum. Inspection of Fig. 1 of Shang et al. (2005) 
suggests that the optical spectrum could have a slightly 
too high flux relative to the HST spectrum. If this is 
indeed the case, then the right-hand-side of 4C+34.47 
extinction curve (beyond 3.5 /im^^) could be lowered by 
2-3 magnitudes and it would then resemble the flatter 
curves in Fig. 4. The reddening, E(B-V), of 4C+34.47 
would then be about 0.03-0.04 mag. less than the values 
in Table 1. The uncertainties in the extinction curves 
for both 4C+34.47 and PG 2349-014 are too great to say 
whether A2175 is present or absent. 

4.3. AGNs With Low Reddening 

For the three AGNs with very low reddening we could 
not get satisfactory extinction curves using all three tem- 
plate curves, so we only used 3C 273. The three extinc- 
tion curves in Fig. 6 are shown for completeness and to 
illustrate that they favor a flatter shape rather than a 
steep SMC-like shape. 

4.4. The Shape of the Far UV Extinction Curve 
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Fig. 6. — Extinction curves for the three AGNs with the lowest 
extinction. UnUke the curves in the two previous figures, these 
curves were derived using only 3C 273 for a template. The smooth 
curves represent SMC, LMC, and Galactic extinction curves. We 
show these curves for completeness and to show that although the 
extinction curves are noisy, (a) they favor a fiatter shape rather 
than an SMC extinction curve shape and (b) there is no rise in the 
far UV. 



One can see in Figs. 4, 5 and 6 that, none of our curves 
show a steep rise towards the far UV. Only PG 0947+396 
shows a shght rise in the far UV. Even the one SMC-hke 
curve, B3 0754+394, seems to flatten in the far UV. Al- 
though we have not derived extinction curves shortwards 
of Lyman a, it is instructive to look at the correlation 
of our reddening values with apuv, the far ultraviolet 
spectral index in a oc v^" power law, meastircd by 
Shang et al. (2005) for A < llOOl. Their auvo, the 
slope from 1200-1500A, is, as would be expected, well 
correlated with our reddenings in Table 1 (correlation 
coefficient, r = 0.87 with a one-tailed probability of aris- 
ing by chance p < 0.001). Since oiuvo is dominated by 
reddening, one might therefore expect apuv to be influ- 
enced by reddening too. As can be seen in Fig. 7, while 
there is a significance correlation (p = 0.04) between the 
reddenings and apuv, f is only 0.45 (i.e., much less than 
for the correlation between auvo and E(B-V).) 

The weak dependence oiapuv on E(B-V) supports the 
idea that, unlike the Galactic curve and the SMC and 
LMC curves, AGN extinction curves do not rise more 
steeply into the far UV. This makes physical sense. The 
rise in far-UV extinction in our Galaxy etc. is the result 
of a broad spectral feature between 10 and 20 eV that 
has a similar origin to that of the A2175 feature (see 
Fig. 2 of Laor & Draine 1993). Laor & Draine (1993) 
consider carbon in the form of graphite, but there is also a 
corresponding spectral feature from PAH molecules (see 
Li & Draine 2001). Since the majority of our extinction 
curves lack obvious A2175, we therefore expect the far 
UV rise to be absent as well. It is interesting the one 
curve showing a rise in the far UV, that of PG 0947+396, 
also shows a possible A2175 feature. 

Although the 10-20 eV peak goes along with the A2175 
bump, the shape of this 10-20 eV feature is much more 
sensitive to the form the carbon is in and the grain size 
distribution than the A2175 feature. Inspection of the 
FUSE composite AGN spectrum of Scott et al. (2004) 




Fig. 7. — The dependence of apjjv, the far ultraviolet spectral 
index, on the reddening for the AGNs of the FUSE/HST sample. 
The line is least-squares fit. 

shows that the sharp peak in the theoretical extinc- 
tion curve of Weingartncr & Draine (2001) due to small 
graphite grains (sec Fig. 9 of Draine 2003) is certainly 
absent. 

The apparently flat far UV extinction curve of AGNs 
does raise one puzzling issue with the far UV: what is 
the cause of the turndown in the spectral shape in the 
far UV? As Shang et al. (2005) show, apuv is steeper 
than auvo by ~ 1. Binette et al. (2005) suggest that 
this is due to the presence of nanodiamonds in AGN 
dust. We suggest, however, that perhaps a more natu- 
ral way of producing this steepening would be by having 
extinction of E(B-V) ~ 0.03 mag caused by intervening 
dust with a Galactic, LMC, or SMC extinction curve, as 
Shang et al. point out. In order not to show up in our 
extinction curves this dust with a steeper far UV extinc- 
tion curve would have to be independent of the nuclear 
dust in the AGN. In agreement with this, our three tem- 
plate AGNs (the rightmost points in Fig. 7) show average 
apuv slopes. 

4.5. A Mean Extinction Curve 

It has long been common practice to deredden AGNs 
with a standard Galactic extinction curve. As we have 
shown above however, AGN reddening curves differ from 
the standard Galactic extinction curve by lacking the 
A2175 feature and by being fiatter in the far UV. In or- 
der to provide a better typical AGN extinction curve for 
de-reddening AGNs we have constructed a mean curve 
from the extinction curves of the AGNs with the high- 
est reddening, but excluding the SMC-like curve of B3 
0754+394. The mean curve was constructed as follows: 
firstly fourth or fifth order polynomials were fit to the 
extinction curves of each of the five AGNs through the 
regions not heavily effected by strong emission lines, and 
then an unweighted average was then taken of these poly- 
nomials. This average extinction curve is shown in Fig. 8 
along with the errors in the mean curve. These errors are 
a minimum around the two regions in which the five in- 
put curves were normalized (see section 3). We also show 
the Galactic extinction curve for the solar neighborhood, 
and the Galactic extinction curve with the Drude pro- 
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Fk;. 8. The average extinction curve for the five AGNs with 
greatest reddening (see text). The mean AGN curve is shown as a 
bold line with 1-sigma error bars for the error in the mean shown 
as dotted lines on either side of it. The solid dashed-dotted curve is 
the Galactic extinction curve for the solar neighborhood, and the 
thin line is the Galactic curve with the A2175 feature removed. 

file of the A2175 feature removed. As can be seen from 

Fig 8, the average AGN extinction curve is in excellent 
agreement with the Galactic curve from the optical up to 
^ 4/im^^ and possibly to 6 or 7 yum"-^, but is significantly 
flatter beyond 7 /um~^ 

For convenience, the mean extinction curve can be rep- 
resented by a fifth-order polynomial m x = in fimr^ , 
valid over the range 1.5 < a; < 8 (i.e., from Ha A6564 to 
Lyman a A 1216). 



£;(A - V)/E{B -V)= 0.000843a;^ - 0.02496a;'^ 
+0.2919a;^ - l.SlSa;^ + 6.83a; - 7.92 

This equation supersedes the relationships given in the 
Appendix of Gaskell et al. (2004). For longer wavelengths 
a standard Galactic extinction curve can be used. 

Although the paramctrization in Eq.(l) gives a more 
appropriate extinction curve than assuming that the 
AGN extinction curve is the same as the Galactic curve, 
researchers arc catitioncd that individual AGN extinc- 
tion curves vary, as can readily be seen in Figs. 4 - 6., 
so one should, if possible, determine the extinction curve 
directly for each individual AGN they are studying. 

5. LUMINOSITY DEPENDENCE 

In Fig. 9 we show the luminosity dependence of our 
reddenings for the FUSE/HST sample compared with 
the luminosity-mean reddening relationship from Gaskell 
et al. (2004). As would be expected, given the bias in 
the FUSE/HST sample towards UV-bright objects, the 
reddenings of the AGNs at a given luminosity mostly 
lie below the typical values found by Gaskell et al. It 
can be seen, however, that the two most reddened AGNs 
for their luminosities have reddenings quite typical of 
the reddenings in the samples in Fig. 5 of Gaskell et al. 
(2004). 

6. CAN REDDENING BE IMITATED BY 
DIFFERENCES IN SPECTRAL-ENERGY 
DISTRIBUTIONS? 




Log(L5ioo) 

Fk;. 9. — The luminosity dependence of the reddenings of the 
AGNs in the FUSE/HST sample. The solid line is from Fig. 5 of 
Gaskell et al. L5100 is XL^ at A5100 in ergs from Shang et al. 
(2005). 

With some notable exceptions (e.g., Carleton et al. 
1987) it has long been widely considered that apparent 
differences in AGN spectral energy distributions (SEDs) 
are caused by intrinsic differences. For example. Yip 
et al. (2004) concluded from their eigenspectral analysis 
of a very large sample of SDSS quasars that the spec- 
tral differences are primarily due to intrinsic differences 
in the AGN spectral energy distribution and to varying 
amounts of host galaxy light. In this section we quanti- 
tatively explore these effects and show that they cannot 
reproduce our observed reddening curves. 

Winkler et al. (1992) find no evidence that the opti- 
cal spectrum changes its shape in an individual AGN as 
it varies, but Barr, Willis & Wilson (1983), and Clavel 
et al. (1991) observe a hardening of the UV spectra of 
some AGNs when they brighten. On the other hand it is 
notable that in the case of Fairall 9 the UV continuum 
slope stayed the same at ajjv = 0.46±0.11 while the UV 
flux changed by a factor of over 20 (Clavel et al. 1989). 
We will model possible spectral shape changes in two 
ways: as a varying "big blue bump" contribution, and as 
a change in the relative contribution of two power-laws. 

6.1. The Effect of Different Big Blue-hump 

Contributions 

Krolik et al. (1991) suggest that the luminosity of the 
"blue bump" varies while an underlying power-law 
remains constant. Such variation in an individual ob- 
ject is inconsistent with the optical and IR observations 
of Winkler et al. (1992), and it is easy to demonstrate 
that such a variation between objects will not explain the 
extinction curves we see. 

In the "power law + blue bump" picture, the blue 
bump is of a thermal origin and has an exponential 
Wien cutoff. We approximated the bump as a single 
black body and added an arbitrary oc i^^^ power law. 
We varied the temperature of the black body, and de- 
rived pseudo-extinction curves from the resulting spec- 
tra. These curves are normalized to E(B-V) = 1 in the 
standard manner . We show our results in Fig. 10 along 
with Milky Way and SMC extinction curves. 

The SED of AGNs probably peaks in the EUV/soft 
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Fig. 10. — Psuedo-extinction curves (thin curves) resulting from 
a change in the temperature of a hot disk compared with Milky 
Way (lower thick curve) and SMC (upper thick curve) extinction 
curves. The top thin curve comes from thermal components of tem- 
peratures 300,000K and 200,000K that are 28% and 45% stronger 
than the power law at Lyman a. The curve just below it is derived 
from components with similar relative strengths and temperatures 
of 100,000K and 50,000K. The dotted curve is derived from com- 
ponents of temperatures 300,000K and 200,000K that are 50 times 
and 80 times stronger than the power law at Lyman a. 



X-ray region of the spectrum (see Gaskell, Klimek & 
Nazarova 2007). This would correspond to tempera- 
tures of 100,000 to 300,000 K. We find that the domi- 
nant parameter affecting the pscudo-rcddcning curves is 
not the temperature or the slope of the arbitrary power 
law, but the ratio of the thermal component to the arbi- 
trary power law. As can be seen from Fig. 10 the choice 
of temperature range has little effect on the shape of the 
pseudo-extinction curves, but varying the ratio of the 
thermal component to the arbitrary power law has a large 
effect. The higher the ratio of the thermal component to 
the arbitrary power law, the flatter the pseudo-extinction 
curve. 

Although it might appear from Fig. 10 that it is easy 
to produce smooth extinction curves resembling those 
we have found, this is very misleading. In reality there 
are two major problems. First, one cannot produce any 
arbitrary reddening. The largest pseudo reddening pro- 
duced in Fig. 10 is only E(B-V) 0.076. Worse still, the 
pseudo reddening is tied to the shape of the extinction 
curve. The second major problem is that the shapes of 
the spectra being compared do not resemble real AGNs. 
In the HST composite of Zheng et al. (1997) the power- 
law index, a, of the spectral region around Lyman a is 
1.72. From the FUSE composite of Scott et al. (2004) a 
= 0.52. For the cases where our pseudo-extinction curves 
resemble the AGN extinction curves derived above, our 
thermal component plus arbitrary power law model gives 
steeply rising spectra with a « -1.7. Such spectra are 
never observed. 

If we try cooler thermal components with temperatures 
of around 30,000 K (see Fig. 11) we can only produce 
SMC- or LMC-like curves. Again, the pseudo-extinction 
curves depend primarily on the strength of the thermal 
component relative to the arbitrary power law and are 
insensitive to the choice of temperatures. It is easier to 
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Fig. 11. — Psuedo-extinction curves (thin curves) resulting from 
a change in the temperature of a cool disk. The thick lines are 
Milky Way and SMC curves as in Fig. 10. Both of the thin 
curves come from comparing thermal components of temperatures 
40,000K and 20,000K. For the upper curve the thermal components 
are 37% and nine times stronger than the power law at Lyman a. 
For the lower thin curve they are 5 times and 37 times stronger 
than the power law at Lyman a. The top and bottom curves cor- 
respond to E(B-V) = 0.098 and 0.175 mag respectively. 

mimic more significant reddenings. However, the spec- 
tral shape problem remains. The 40,000K black body 
plus power law produces gives a rising spectrum of a ~ - 
0.5 in the spectral region around Lyman a. This is again 
at variance with observations. 

6.2. Two Power Laws? 

We also considered an arbitrary mixture of power-laws. 
The overall optical to X-ray spectral shape of most AGNs 
can be described with an a = 1.7 power law. Winkler et 
al. (1992) present the shapes of the variable continua for 
many AGNs. We took the bluest of these to represent 
the unreddened variable continuum and represent this 
with an a = -0.2 power law. In Fig. 12 we show the 
pseudo-extinction curves obtained by adding increasing 
amounts of the a = -0.2 power law to the a = 1.7 power 
law. The precise choice of spectral indices of the two 
components is not important. As can be seen from Fig. 
12, it is possible to reproduce SMC- or LMC-like curves 
and somewhat higher apparent reddenings, but the same 
problems exist as with the thermal components in Fig. 
11: the shape of the extinction curve is coupled to the 
degree of reddening (see caption to Fig. 11), it is not 
possible to get very high reddenings, and for all but the 
most modest reddenings, the shapes of the bluest spectra 
are again much steeper than is observed. In Fig. 13 we 
show the continuum spectral indices around Lyman a for 
the bluest spectra as a function of the pseudo reddening. 

6.3. The Effect of Differing Host Galaxy Contributions 

The final effect we consider is the idea that the differ- 
ences in spectral energy distributions are due to chang- 
ing host galaxy contributions (Yip et al. 2004). We used 
a simple cool black body to approximate the effect of 
a changing host galaxy stellar population and approxi- 
mated the quasar continuum as a power law. As 
is to be expected, all the changes are in the long wave- 
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Fig. 12. — Psuedo-extinction curves (thin curves) resulting from 
adding varying amounts of an a = -0.2 power law to an a = 1.7 
power law. The thick lines are Milky Way and SMC curves as in 
Figs. 10 and 11. The thin lines correspond to effective reddenings 
of E(B-V) = 0.02, 0.05, 0.13, 0.2.5, 0.36, & 0.43 magnitudes. 
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Fig. 13. — Continuum spectral index in the Lyman a spectral 
region for the bluest AGNs as a function of the pseudo reddening, 

E(B-V)e//. 

length end of the spectrum, and if the changes in spectral 
shape are interpreted as a pseudo-extinction curve (see 
Fig. 14), the UV part of the curve bears no resemblance 
to a reddening curve at all. 

From the analysis presented in this section we con- 
clude that changes in the continuum shape from AGN 
to AGN, either because of changes in the AGN itself or 
in the degree of host galaxy contamination, are not the 
cause of the extinction curves we have derived. This con- 
clusion is also supported by the similarity of broad-line 
and continuum reddenings, and reddenings derived from 
X-rays (see discussion in section 7.6). A positive conclu- 
sion that can be drawn from the analysis of sections 6.1 
and 6.2, however, is that because changes in continuum 
shape do mimic reddening to some degree, small changes 
in the shape of the continua from object to object, if 
such changes exist, will not be significantly distorting 
the shapes we derive. 

7. DISCUSSION 

7.1. The Fraction of Steep Curves 

Richards et al. (2003) only consider the possibility of 
AGNs being reddening by SMG-like dust, and Willott 
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Fig. 14. — Similar to Figs. 11 and 12, but showing the effect of 
a cooler black body changing its temperature from 4000 to 5000 
K (i.e., the range of effective temperature of an old stellar popula- 
tion), and its contribution to the flux at 5500A from 19% to 72% 
compared with a a = — 1 power law. 

(2005) also argues for AGNs being reddened by SMC- 
like dust, but in our study we find only one clear case of 
a steep, SMC-hke reddening curve, B3 0754+394. We be- 
lieve the paucity of SMC-likc curves is supported by other 
studies. Gaskell et al. (2004) found quite flat extinction 
in the UV for the samples of AGNs they considered, espe- 
cially the radio-loud AGNs. Crenshaw et al. (2001) found 
a clearly steep, SMC-likc extinction curve for the radio- 
quiet Seyfert NGC 3227 (E(B-V) « 0.180 compared with 
NGC 4151). However, the extinction curve Crenshaw et 
al. (2002) found for Ark 564 is similar to a number of 
those in Figs. 4 and 5. Gaskell, Klimek, & Nazarova 
(2007) find an extinction curve for NGC 5548 similar to 
our average curve presented here. 

Richards et al. (2003) only consider SMC-like redden- 
ing curves. Gaskell et al. (2004) have pointed out that the 
statement by Richards et al. that attempting to explain 
their large color-segregated composites with reddening 
"results in good matches at both 1700 and 4040 A but 
overpredicts the flux between these two wavelengths and 
underpredicts the flux shortward of C IV" shows that the 
reddening curve for the AGNs they consider must in fact 
be flatter than an SMC-like curve. 

As noted in section 2, the FUSE/HST/optical sam- 
ple of Shang et al. (2005) is biased towards AGNs with 
low reddening. More specifically, it is biased towards 
AGNs with low ultraviolet reddening. It will therefore be 
particularly biased against AGNs which are reddened by 
SMC-like extinction curves. We therefore predict that a 
less biased sample would not only have more AGNs with 
higher reddening, but could also have more with steeper 
(e.g., SMC-like) extinction curves. However, the small 
spread in the UV slope of AGNs (Cheng, Gaskell, & Ko- 
ratkar 1991) suggests that steep, SMC-like curves are not 
common. 

7.2. The X2175 Feature 

The A2175 feature has long been known to be weak 
in AGNs (McKee & Petrosian 1974). Of the six highest 
quality extinction curves presented here (see Fig. 4), none 
show the A2175 feature, while there is an indication that 
three of the next highest quality five extinction curves 
might show it. Nonetheless we believe that the A2175 fea- 



AGN EXTINCTION CURVES 



9 



ture is relatively rare in AGN extinction curves because 
it is completely absent in the extinction curves presented 
by Czerny et al. (2004) and Gaskell et al. (2004) based on 
composite spectra. Although, as noted in the introduc- 
tion, luminosity biases can make extinction curves from 
composite spectra (e.g., those of Gaskell et al. 2004, and 
Czerny ct al. 2004) flatter in the UV than they ought to 
be (Willott 2005), one thing this will not do is affect the 
strength of A2175. 

7.3. Emission-line Reddenings 

Emission-line reddenings can be estimated from known 
line ratios. The He II A1640/A4686 and O I A1304/A8446 
ratios are particularly useful. For Mrk 290 the redshift 
is low enough that O I A8446 was observed, and the op- 
tical Fe II emission is weak enough that He II A4686 can 
be measured. Relative line strengths were estimated by 
profile fitting. We fit the He II lines with the C IV profile 
and the O I lines with the Ha profile. We estimate that 
the line ratios are uncertain by up to about 50%. This is 
because the lines are weak and determining the appropri- 
ate underlying continuum in the presence of other weak 
lines and the wings of stronger lines is difhcult. For Mrk 
290 we estimate the He II A1640/A4686 ratio to be 3.8. 
If the real ratio is 8.3, as suggested by photoionization 
models, then our average extinction curve implies a total 
reddening, E(B-V) of 0.20, or an intrinsic reddening of 
0.18 after allowance for Galactic reddening of 0.016 mag. 
We estimate the O I AI304/A8446 ratio to be 3.1 com- 
pared with the theoretical value of 6.5. This implies an 
intrinsic reddening of 0.17 mag. 

These two reddening estimates are in agreement with 
our estimate of 0.16 mag from the continuum shape (i.e., 
the; value in Table 1). The errors in the reddenings es- 
timated from the O I and He II ratios could be up to 
±0.1 mag, but the agreement provides support for the 
differences in observed continuum shapes being due to 
reddening. Further careful study of these line ratios in 
more AGNs would be very worthwhile to check what the 
zero-point of the reddening is (i.e., what the reddening 
is for our template AGNs). 

7.4. The Shape of the Spectral- Energy Distribution 

Emission lines provide additional support for our con- 
tention that apparent differences in the optical-UV spec- 
tral shape are first and foremost due to reddening. We 
have shown that reddening, and hence the apparent op- 
tical to UV spectral index, is luminosity dependent (see 
Gaskell et al. 2004 and section 9 above). There is no 
doubt that broad emission lines are produced by pho- 
toionization resulting from an ionizing continuum that 
is an extrapolation of the observed UV-optical contin- 
uum. If the UV-optical continuum really were steeper 
at lower luminosities, rather than being reddened more, 
then there would be far fewer ionizing photons on aver- 
age in lower-luminosity AGNs. This would result in much 
lower emission-line equivalent widths in lower-luminosity 
AGNs, which is completely the opposite of the Baldwin 
effect. 

With the possible exception of the far UV (see sec- 
tion 4.4), in the spectral region we consider, the only 
spectral differences we detect that are clearly not due to 
reddening are due to differing broad-line region proper- 
ties. These include differing strengths and widths of the 



major emission lines, and differences in the "small blue 
bump" due to Balmer continuum and Fe II emission. 

7.5. Implications of Flat Versus Steep Reddening 
Curves 

The shape of the reddening curve strongly affects the 
estimated extinction. For a given extinction at A1640, 
say, a flat curve such as that in Fig. 8 implies an E(B- 
V) about three times greater than what an SMC curve 
would imply. This means that by assuming an SMC-like 
curve a priori, Richards et al. (2003) are underestimating 
the reddening of SDSS AGNs (e.g., in their Fig. 6) by a 
factor of 3-5 at higher redshifts. 

The larger grain sizes producing the flatter cure could 
also mean that the ratio, R, of V-band extinction, Ay, to 
E(B-V) is greater than the value of 3.1 commonly found 
in the local ISM. Gaskell et al. (2004) suggest that i? « 5 
is appropriate for their flat extinction curve. 

Taken together, the increased E(B-V) and, to a lesser 
extent, the possibility of an increased R, means that Ay 
is a lot greater with a flatter extinction curve than with 
an SMC-like curve. This has important implications for 
AGN energetics and demographics. 

7.6. The Amount of Reddening of AGNs 

De Zotti & Gaskell (1985) flnd, from considerations 
of broad-band optical colors and Balmer decrements as 
a function of orientation, that the typical reddening for 
a large sample of optically-selected Seyfert 1 galaxies is 
~ 0.30 mag. Winkler et al. (1992) derive reddenings for 
a large sample of southern Scyfcrts from Balmer decre- 
ments, the X-ray to H/3 ratio, and variability colors. The 
mean reddenings from their three methods are E(B-V) = 
0.32, 0.41, and 0.29. From a variety of broad-line ratios 
Ward & Morris (1984) get E(B-V) = 0.27 mag. for NGC 
3783. Estimates of the reddening from the UV/optical 
continuum shapes of individual AGNs are consistent with 
these reddenings. For NGC 3227 Crenshaw et al. (2001) 
get E(B-V) = 0.18 relative to NGC 4151, and for Ark 
564 Crenshaw et al. (2002) get E(B-V) = 0.14 compared 
with Mrk 493. For NGC 5548 Gaskell et al. (2005) get 
E(B-V) = 0.20 compared with 3C 273. AU of these AGNs 
were optically selected. Considering all of these redden- 
ing estimates we see that, as expected, the reddenings of 
the FUSE/HST/optical sample are lower than average 
for AGNs. 

When one looks at complete samples of AGNs se- 
lected by radio emission, the reddenings are even greater. 
For the four increasingly reddened Molonglo samples of 
Baker & Hunstead (1995) Gaskell et al. (2004) get mean 
reddenings of E(B-V) = 0.29, 0.34, 0.71, and 0.98 mag. 
Thus the Baker & Hunstead samples are reddened more 
than typical optically-selected Seyfert galaxies, and the 
least reddened of the Molonglo samples is more red- 
dened than the FUSE/HST/optical sample of Shang et 
al. (2005). 

From all these diverse reddening estimates it is diffi- 
cult to avoid the conclusion that reddening must have a 
strong effect on the apparent spectral shapes of AGNs. 
Additional, independent support for steeper slopes being 
a consequence of dust comes from the correlation between 
spectral index and the strength of associate absorption 
(Wills, Shang, & Yuan 2000, Baker et al. 2002). 
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8. CONCLUSIONS 

In agreement with Gaskell et al. (2004) we argue that 
the UV-optical continuum shape is profoundly affected 
by reddening for all but the bluest AGNs. We believe 
that all apparent differences in the UV-optical continua 
of AGNs are primarily due to differences in extinction 
and in the forms of the extinction curves. Until the 
effects of reddening are carefully allowed for it is not 
possible to investigate possible remaining differences in 
the spectral energy distributions of AGNs as a function 
of such things as luminosity, black-hole mass, accretion 
rate, or radio power. At present all AGNs have optical- 
UV continuous spectra that appear to be consistent with 
a single shape after allowance for the effects of reddening 
and the small blue bump. Even in a sample, such at the 
HST/FUSE sample, biased towards low reddening, most 
of the objects show detectable reddening. 

Our results support the interpretation of Shang et al. 
(2003) that the second eigenvector in their spectral prin- 
cipal component analysis is due to dust reddening, and 
disagree with the interpretation of Yip et al. (2004) that 
the main continuum eigenvectors are due to changing 
host galaxy contributions and variations in the intrinsic 
spectral index. 

We find rising SMC-like reddening curves to be rare 
in the relatively low reddening HST/FUSE sample we 
consider and believe that they are not the most common 
reddening curve for AGNs. Flatter reddening curves im- 
ply that the visual extinction of AGNs is greater than if 
the reddening curves were SMC-like. 



We confirm that the A2175 extinction feature is mostly 
absent in AGN extinction curves. The rise into the far 
UV seen in the Galactic, LMC, and SMC reddening 
curves seems also to be mostly absent. These differences 
imply that the dust has been modified by the AGN en- 
vironment. 

Our new mean AGN extinction curve, while flat in 
the far UV (1/A > 6.5 /im~^) is steeper in the range 
3.5 < 1/A < 6.5 /xm-i than the Gaskell et al. (2004) 
extinction curve from composite spectra of radio-loud 
AGNs. Although the objects are not the same, the 
difference suggests that himinosity-dependent extinction 
(Willott 2005) might have caused the Gaskell et al. curve 
to be too flat in the UV. In disagreement with Willott 
(2005), however, we do not find the majority of AGN 
extinction curves to be SMC-like. 
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